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We present the next-to-leading order (NLO) QCD results for W + 4-jet production at hadron 
colliders. This is the first hadron collider process with five final-state objects to be computed 
at NLO. It represents an important background to many searches for new physics at the energy 
frontier. Total cross sections, as well as distributions in the jet transverse momenta and in the total 
transverse energy Ht, are provided for the initial LHC energy of y/s = 7 TeV. We use a leading-color 
approximation, known to be accurate to 3% for W production with fewer jets. The calculation uses 
the BlackHat library along with the SHERPA package. 
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The first data and analyses emerging from experiments 
at the Large Hadron Collider (LHC) emphasize the need 
for reliable theoretical calculations in searches for new 
physics beyond the Standard Model. In many channels, 
new-physics signals can hide in broad distributions un- 
derneath Standard Model backgrounds. Extraction of 
a signal will require accurate predictions for the back- 
ground processes, for which next-to-leading order (NLO) 
cross sections in perturbative QCD are crucial. The past 
few years have seen rapid progress in NLO QCD for the 
LHC. Several important processes involving four final- 
state objects (including jets) have been computed [H-EI . 

In this Letter, we present results for the first of a new 
class of processes, involving five final-state objects: inclu- 
sive W + 4-jet production, using a leading-color approx- 
imation for the virtual terms that has been validated for 
processes with fewer jets. This process has been studied 
since the early days of the Tevatron, where it was the 
dominant background to top-quark pair production. At 
the LHC, it will be an important background to many 
new physics searches involving missing energy, as well as 
to precise top-quark measurements. 

In previous papers [H [2| we presented the first com- 
plete results for hadron-collidcr production of a W or 
Z boson in association with three jets at NLO in the 
strong coupling a s . (Other NLO results for W + 3 jets 
have used various leading-color approximations [2Hs}L) 
We performed detailed comparisons to Tevatron data [7] . 





FIG. 1: Sample diagrams for the seven-point loop amplitudes 
for qg — » W q'ggg and qQ — > W q'ggQ, followed by W — > ev. 



The sensitivity to the unphysical scale used to define a s 
and the parton distributions is reduced from around 40% 
at leading order (LO) to 10-20% at NLO, and the NLO 
results agree well with the data. At the LHC, a much 
wider range of kinematics will be probed, making NLO 
studies even more important. 

The computation of hadron collider processes with 
complex final states at NLO has long been a challenge to 
theorists. The evaluation of the one-loop (virtual) cor- 
rections has been a longstanding bottleneck. Feynman- 
diagram techniques suffer from rapid growth in complex- 
ity as the number of legs increases. On-shell methods [§- 
[l3j , in contrast, rely on the unitarity and factorization 
properties of scattering amplitudes to generate new am- 
plitudes from previously-computed ones. Such methods 
scale very well as the number of external legs increases, 
offering a solution to these difficulties. 
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We use the same basic setup as in our earlier computa- 
tions of W+ 3-jet and Z, 7* + 3-jet production. The 
virtual contributions are computed using on-shell meth- 
ods via the BlackHat package [14j . We show represen- 
tative virtual diagrams in fig. [TJ We use a leading-color 
approximation in the finite virtual contributions 1 , while 
keeping the full color dependence in all other contribu- 
tions. We have confirmed that this approximation is an 
excellent one for W + 1,2, 3-jet production, shifting the 
total cross section by about 3%, which is significantly 
smaller than uncertainties from parton distributions or 
higher-order terms in a s . Subleading-color corrections to 
W + 4-jet production should also be small. We include 
the full W Breit-Wigner resonance; decays to leptons re- 
tain all spin correlations. 

The remaining NLO ingredients, the real-emission 
and dipole-subtraction terms [Hj], are computed by 
AMEGIC++ 0, part of the SHERPA package (nf. 
We also use SHERPA to perform phase-space integra- 
tion. The efficiency of the integrator has been improved 
significantly with respect to ref. [l| through the use of 
QCD antenna structures 3, 19|. BlackHat computes 
the real-emission tree amplitudes using on-shell recur- 
sion relations (lo| , along with efficient analytic forms ex- 
tracted from M = 4 super- Yang-Mills theory [2(| ■ 

Compared to LO, NLO cross sections and distributions 
generally depend much less on the common (unphysical) 
renormalization and factorization scale /x. However, it is 
still important to select a scale characteristic of the typ- 
ical kinematics. A scale that performs well for many dis- 
tributions is the total partonic transverse energy. We set 
H = H' T I '2, where H' T = Y]- pip + E^f ; the sum runs over 
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(pf) 2 is the 



all final-state partons j, and E, 

transverse energy of the W boson 2 . Refs. [1, [22| present 
other satisfactory choices. We follow the conventional 
procedure of varying the chosen central scale up and 
down by a factor of two to construct scale-dependence 
bands, taking the minimum and maximum of the ob- 
servable evaluated at five values: fi/2, (J,, \2~li, 2/i. 

Fixed-order perturbation theory may break down in 
special kinematic regions, where large logarithms of scale 
ratios emerge. For instance, threshold logarithms can af- 
fect production at very large mass scales, which can be 
reached in inclusive single-jet production 23[ . Using this 
study one can argue [2| that at the mass scales probed 
in W + 4-jet production, such logarithms should remain 
quite modest. Similarly, the sort of large logarithms aris- 
ing in vector-boson production in association with a sin- 



1 Our definition of leading-color terms follows ref. Q|; it includes 
virtual quark loops in addition to the terms identified in ref. 

2 In refs. Q, we used the scalar sum of the decay leptons' trans- 
verse energies instead of . The present choice is preferred for 



gle jet [24[ do not appear in the case of multiple jets. 
Tighter cuts may isolate regions which would require a 
reassessment of potentially-large logarithms. 

In our study, we consider the inclusive process pp — > 
W + 4 jets at an LHC center-of-mass energy of </s = 
7 TeV. We impose the following cuts: E^ > 



20 

,jet| 



studies of W polarization effects [1|, 12. 



GeV, 

|?f I < 2.5, $ T > 20 GeV, p 3 ? > 25 GeV, |^ Gt | < 3, 
and Mj[ > 20 GeV. Here, pr are transverse momenta; 
77, pseudorapidities; and My, the transverse mass of the 
ev pair. The missing transverse energy, corresponds 
to the neutrino transverse energy, Ey. Jets are defined 
using the anti-fey algorithm [25| with parameter R = 
0.5, and are ordered in py. (We also quote results for 
R = 0.4.) We use the CTEQ6M [H parton distribution 
functions and a s at NLO, and the CTEQ6L1 set at LO. 
Electroweak couplings are as in ref. [l[. 

In table UJ we present LO and NLO parton- level cross 
sections for inclusive W^-boson production accompanied 
by zero through four jets. We include all subprocesses, 
using the leading-color virtual approximation only in W+ 
4-jet production. The upward scale-variation figures for 
the NLO cross sections are quite small for W + 3- and 
W + 4-jet production, because the values at the central 
scale choice are close to the maximum values across scale 
variations. We also display the ratios of the W + to W~ 
cross sections, and the "jet-production" ratios of W~ + n- 
jet to W~ + (n— l)-jet production. Both kinds of ratios 
should be less sensitive to experimental and theoretical 
systematics than the absolute cross sections. 

The W + fW~ ratios are greater than one because the 
LHC is a pp machine, and because the parton luminos- 
ity ratio u(x)/d(x) exceeds one. As the number of jets 
increases, production of a W requires a larger value of x, 
driving u(x)/d(x) and hence the W + jW~ ratio to larger 
values. These ratios have been discussed recently 27j as 
a probe of certain new-physics processes; our results ex- 
tend the NLO analysis to W production accompanied by 
four jets. This ratio changes very little under correlated 
variations of the scale in numerator and denominator; 
hence we do not exhibit such scale variation here. 

Standard lore [28j says that the jet-production ratio 
should be roughly independent of the number of jets. The 
results for the ratios displayed here for n > 1 are indeed 
consistent with this lore. However, they are rather sensi- 
tive to the experimental cuts, and can depend strongly on 
n when binned in the vector-boson px [2| . The W + 1- 
jet/iy + 0-jet ratio is much smaller because of the re- 
stricted kinematics of the leading contribution to W + 0- 
jet production. 

In table [Til we g ive cross sections for narrower jets, 
using the anti-fey jet algorithm with R = 0.4. For two 
or more jets, the LO cross sections are larger than for 
R = 0.5, and the effect increases with the number of jets. 
However, at NLO, the effect is greatly diminished; only 
for four jets is the NLO cross section for R — 0.4 signifi- 
cantly above that for R = 0.5. The NLO jet-production 
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FIG. 2: A comparison of the pr distributions of the leading four jets in W~ + 4-jet production at the LHC at \fs = 7 TeV. 
In the upper panels the NLO distribution is the solid (black) histogram and the LO predictions are shown as dashed (blue) 
lines. The thin vertical line in the center of each bin (where visible) gives its numerical (Monte Carlo) integration error. The 
lower panels show the LO distribution and LO and NLO scale-dependence bands normalized to the central NLO prediction. 
The bands are shaded (gray) for NLO and cross-hatched (brown) for LO. 



no. jets 


W~ LO 


W~ NLO 


W+/W~ LO 


W+/W~ NLO 


W-n/(n-l) LO 


W~n/(n-l) NLO 





1614.0(0.5)1^ 


2077(2) +*? 


1.656(0.001) 


1.580(0.004) 






1 


264.4(0.2)±^;| 


331(1)13 


1.507(0.002) 


1.498(0.009) 


0.1638(0.0001)1^31 


0.159(0.001) 


2 


73.14(0.09)+^ 


78.1(0.5)1^ 


1.596(0.003) 


1.57(0.02) 


0.2766(0.0004)1°;°^ 


0.236(0.002) 


3 


17.22(0.03)1^ 


16.9(0.1)1°J 


1.694(0.005) 


1.66(0.02) 


0.2354(0.0005)l°;o25 


0.216(0.002) 


4 


3.81(0.01)1?;** 


3.55(0.04)1^ 


1.812(0.001) 


1.73(0.03) 


0.2212(0.0004)1^20 


0.210(0.003) 



TABLE I: Total cross sections in pb for W + n jet production at the LHC at ^fs = 7 TeV, using the anti-fcr jet algorithm 
with R — 0.5. The NLO result for W + 4 jets uses the leading-color approximation discussed in the text. The fourth and 
fifth columns give the cross-section ratios for W + production to W~ production. The last column gives the ratios of the cross 
section for the given process to that with one jet less. The numerical integration uncertainty is in parentheses, and the scale 
dependence is quoted in super- and subscripts. 



ratio is somewhat larger for R = 0.4, for n > 2; in con- 
trast, the ratios of W + to W~ cross sections are un- 
changed within errors. 

In fig. [21 we show the px distributions of the leading 
four jets in W~ + 4-jet production at LO and NLO; the 
predictions are normalized to the central NLO prediction 
in the lower panels. With our central scale choice, there is 
a noticeable shape difference between the LO and NLO 
distributions for the first three leading jets, while the 
fourth-jet distribution is very similar at LO and NLO. 
Similarly, in W + 3-jet production, the pr distributions 
of the leading two jets exhibit shape changes from LO to 
NLO, while the third-jet distribution does not [l|. 

Fig. [3] shows the distribution of the total transverse 
energy Ht, given by the scalar sum of the jet and lep- 
ton transverse energies, Ht = J2j -^t*" + + ^T- We 



show the NLO and LO predictions, along with their scale- 
dependence bands. As in the px distributions, the NLO 
band is narrower. The shapes at LO and NLO are similar 
above 200 GeV, where the integration errors are small. 

The results of this study validate our understanding of 
the W + 4-jet process for typical Standard-Model cuts. 
It will be interesting, and necessary, to explore the size of 
corrections for observables and cuts used in new-physics 
searches. 

In order to compare our parton-level results to forth- 
coming experimental data, the size of non-perturbative 
effects (such as hadronization and the underlying event) 
needs to be estimated, for example using LO parton- 
shower Monte Carlo programs. As NLO parton-shower 
programs are developed [29| , the virtual corrections com- 
puted here should be incorporated into them. 
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no. jets 


W~ LO 


W~ NLO 


W + /W~ LO 


W + /W~ NLO 


W-n/(n-l) LO 


W~n/(n-l) NLO 





1614.0(0.5)™ 


2077(2)1* 


1.656(0.001) 


1.580(0.004) 






1 


264.4(0.2)±^ 


324(1)1}? 


1.507(0.002) 


1.499(0.009) 


0.1638(0.0001)1°,;°,^ 


0.156(0.001) 


2 


74.17(0.09)1^ 


76.2(0.5)13.4 


1.597(0.003) 


1.56(0.02) 


0.2805 (0.0004) looas 


0.235(0.002) 


3 


18.42(0.03)lg;2g 


17.0(0.1)1?° 


1.694(0.005) 


1.66(0.02) 


0.2483(0.0005)1q;o26 


0.223(0.002) 


4 


4.41(0.01)1?;^ 


3.81(0.04)1°°° 


1.814(0.001) 


1.76(0.03) 


0.2394(0.0004)1^2? 


0.224(0.003) 



TABLE II: The same quantities as in table HI but with R = 0.4. 



200 300 400 500 600 700 800 900 1000 




LO / NLO I NLO scale dependence LO scale dependence 



1.5 L 

: |pf ' ' 

0.5 R"^ i .... i i .... i ,,, , — 

200 300 400 500 600 700 800 900 1000 
H T [GeV] 

FIG. 3: The H T distribution for W~ + 4 jets. 



A related process that contributes an irreducible back- 
ground to certain missing energy signals of new physics is 
Z+ 4-jet production. We expect that the current Black- 
Hat along with SHERPA will allow us to compute NLO 
corrections to it, as well as to other complex processes, 
thereby providing an unprecedented level of theoretical 
precision for such backgrounds at the LHC. 
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